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Abstract. We recover spiral and barred spiral patterns in disk galaxy simnlations with a 
Wave Dark Matter (WDM) backgronnd (also known as Scalar Field Dark Matter (SFDM), 
Ultra-Light Axion (ULA) dark matter, and Bose-Einstein Condensate (BEC) dark matter). 
Here we show how the interaction between a baryonic disk and its Dark Matter Halo triggers 
the formation of spiral strnctnres when the halo is allowed to have a triaxial shape and angn- 
lar momentnm. This is a more realistic pictnre within the WDM model since a non-spherical 
rotating halo seems to be more natnral. By performing hydrodynamic simnlations, along with 
earlier test particles simnlations, we demonstrate another important way in which wave dark 
matter is consistent with observations. The common existence of bars in these simnlations 
is particnlarly noteworthy. This may have conseqnences when trying to obtain information 
abont the dark matter distribntion in a galaxy, the mere presence of spiral arms or a bar nsn- 
ally indicates that baryonic matter dominates the central region and therefore observations, 
like rotation cnrves, may not tell ns what the DM distribntion is at the halo center. Bnt 
here we show that spiral arms and bars can develop in DM dominated galaxies with a central 
density core withont snpposing its origin on mechanisms intrinsic to the baryonic matter. 
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1 Introduction 

There are many explanations on the origin of the Spiral structure observed in disk galaxies. 
It is accepted that some of these spiral structures are driven by tides [1], by a response to 
bars [2], or induced by density clumps such as giant molecular clouds [3]. But it is also a 
result of numerical simulations that spiral patterns can be self-excited structures, developing 
without the need of a trigger. In general the spiral arms resulting from N-body simulations 
are multiple, transient features [4-6] that span a range of pattern speeds, but there is no 
agreement on the lifetimes of the patterns. 

Aside from the trigger mechanisms mentioned above, or the complexity of self-excited 
spirals seen in simulations, the interaction between a baryonic disk and its Dark Matter Halo 
has been studied for a few decades. It has been proposed to be responsible for a warp formation 
[7, 8]. A non-spherical halo misaligned with the disk imprints a torque that could form a warp. 
Because dark matter halos in cosmological simulation are non-spherical, with misalignments 
between the principal axes of the inner and outer halo, Dubinski & Chakrabarty [9] obtain 
warps in their simulations when the disk is perturbed by such slowly rotating, rigid, halos. 
Following this approach, it is plausible to consider a triaxial dark matter halo as the trigger 
mechanism of spiral structure formation in the disk. Halos with strongly triaxial shapes are 
predicted by cosmological simulations within the Cold Dark Matter (CDM) scenario [10, 11] 
via the formation process of galactic size halos occurring through the hierarchical assembly 
of CDM structures [12, 13]. 

The formation of a spiral pattern in the presence of a triaxial dark matter halo has been 
studied in several works. For example the spiral structure extended beyond optical limits in 
the dwarf galaxy NGC 2915 can be explained by the presence of a triaxial dark halo [14] A 
variety of halo density profiles have been used in these works, such as the NFW profile [14], a 
logarithmic halo [15], or a quasi-isothermal halo [16]. The common result of these studies is 
that, in spite the different dark halo models, formation of spiral structure in the stellar and 
gaseous disk is induced by the presence of a dark halo with triaxial shape. 

The works mentioned above, that use N-body or hydrodynamics to simulate the stellar 
and gas components, model the DM halo as a rigid triaxial potential. By starting with a 
spherically symmetric DM distribution, the usual approach to account for triaxiality is to 
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generalize the isodensity snrfaces into triaxial ellipsoids. Here we also stndy the interaction 
between the DM halo and the disk of the galaxy by nsing a rigid halo, bnt acconnt for 
triaxiality not following the nsnal approach, bnt in a more natnral way. This paper focnses 
on Wave Dark Matter which we snggest as an nmbrella term for Scalar Field Dark Matter 
(SFDM), Ultra-Light Axion (ULA) dark matter, and Bose-Einstein Condensate (BEC) dark 
matter becanse of the wave natnre of these nearly identical models on the scale of galaxies. 

Wave dark matter and the CDM model make virtnally identical predictions for the 
distribntion of dark matter at cosmological scales ([17] or see recent nnmerical simnlations 
confirming this resnlt in [18]). Bnt is worth mention that the cnspy density profiles of halos 
in CDM-only simnlations snppose a contradiction with high qnality observations, a core-like 
behavior is preferred in the central regions of dwarfs Spheroidals (dSphs) [19, 20] and Low 
snrface brightness galaxies (LSB) [21-23]. CDM faces these and other difhcnlties on galactic 
scales, it appears to be necessary to look at the predictions of alternative dark matter models 
at scales of galaxies to overcome these difhcnlties. In onr alternative scenario, dark matter 
is modeled by a scalar held satisfying the Klein-Gordon eqnation [24-28], a wave type of 
eqnation. These types of wave eqnations, throngh constructive and destrnctive interference, 
natnrally form density waves which, throngh gravity, may indnce density waves in the visible 
regnlar matter of a galaxy. Onr simnlations demonstrate that these visible density waves 
often form spiral and barred spiral patterns very similar to those observed in the nniverse. 
Hence, the fact that wave dark matter easily prodnces spiral and barred spiral patterns in 
galaxies is very relevant to the problem of determining the model that better captnres the 
natnre of dark matter. 

The ontcome that a scalar held conhgnration has a spacial distribntion that in general 
is triaxial, triggering the formation of spiral and barred patterns, reinforce the hypothesis of 
a scalar held as dark matter. This reveals another important galaxy-like behavior of a scalar 
held conhgnration, along with the ability to model galactic systems (Lsb and dSph galaxies), 
as tested with observational data [29-31]. 

We comment that the recent discovery of the Higgs boson which is modeled by a scalar 
held as well, snpports the idea that scalar helds may indeed be physical. One interpretation 
of the wave dark matter model is that dark matter is yet another spin zero boson, with a mass 
on the order of 10“^^ eV to 10“^*^ eV [17, 32, 33] and with wavelengths perhaps approaching 
the scale of the smallest galaxies. In this regime, the wave-particle dnality of matter canses 
the boson to act like a classical scalar held satisfying a wave type of eqnation, namely the 
Klein-Gordon eqnation. 

This paper is organized as follows. First, in Sect. 2 we present the wave dark matter 
triaxial potential. Then, in order to compare with another triaxial potential, in Sect. 3 we 
present a Navarro-Frenk-White triaxial halo. In Sect. 4 we briehy describe the nnmerical code 
and initial conditions for onr simnlations. In Sect. 5 we present the resnlts of the simnlations 
with the hydrodynamic code and triaxial halos. Finally, a discnssion and conclnsions are 
presented in Sect. 6. 

2 Wave Dark Matter in a Spiral Galaxy 

We model the wave dark matter in this paper precisely as in Bray (2010). We refer the 
reader to that paper for details and give an overview in this paper. To hnd approximate 
solntions to the nonlinear Einstein-Klein-Gordon eqnations, we consider exact solntions to 
the linear Klein-Gordon eqnation in a Minkowski spacetime backgronnd. We then cnt these 
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linear solutions off at an arbitrary radius. The effect of this is to produce what we claim are 
qualitatively realistic solutions to the nonlinear Einstein-Klein-Gordon equations, as explained 
in Bray (2010). 

The linear solution we choose is the sum of a spherically symmetric term and a second 
degree spherical harmonic term. These two terms interfere, producing an interesting wave 
pattern which rotates around the center of the galaxy. Since we’re in the Newtonian regime, 
we may use the Poisson equation to recover the corresponding Newtonian potential which 
is also rotating and is ellipsoidal. It is precisely this rotating ellipsoidal galactic potential - 
qualitatively predicted by the wave dark matter model - that we feed into the simulations of 
the stars, gas, and dust. 

The addition of the second degree spherical harmonic term to the scalar held allows the 
wave dark matter to have angular momentum, for example, which seems natural. We do not 
add in a hrst degree spherical harmonic term because this would cause the center of mass of 
the wave dark matter, and hence the regular matter, to oscillate. We believe the friction of the 
gas and dust of spiral galaxies would damp out this behavior. We do not add in higher degree 
terms because they drop off more quickly than the second degree term and are assumed to be 
less relevant in most situations. Future simulations could add in more terms to produce more 
interesting features. In particular, galaxies with three or more spiral arms could be related 
to instances when third degree or higher terms dominate the second degree term. 


2.1 Triaxial Potential 

As explained in Bray (2010), wave dark matter and its associated gravitational held are 
approximately modeled by the Poisson-Klein-Gordon equations: 






f 




( 2 . 1 ) 

( 2 . 2 ) 


where V is the Newtonian gravitational potential of the galaxy, fioM is the wave dark mat¬ 
ter density, /tq is the density of the wave dark matter at the origin times the Newtonian 
gravitational constant G, ^ is the standard Laplacian, f{t, x, y, z) is a 

function (or scalar held) representing the dark matter, and m is the mass of the dark matter 
boson being modeled times c/h. Equation 2.2 is a good approximation for long wavelength 
solutions for /, meaning that / is changing in time much more quickly than in space. The 
main physical characteristic that we have lost in this system is that V(t, x, y, z) does not feed 
back to affect f{t,x,y,z) in the second equation. In other words, it appears that we are not 
accounting for the gravity of the galaxy on the scalar held. Of course, gravity is what keeps 
the scalar held (and everything else) trapped in the galaxy. But we can qualitatively account 
for this missing effect by simply requiring that the scalar held / goes to zero outside of some 
large radius. 

In other words, we are accounting for the effect of gravity (from all sources) on the 
scalar held / representing dark matter, at least up to a fairly good approximation. If the 
total gravitational held is assumed to be spherically symmetric in the galaxy (which it roughly 
is), then one can study the effect this gravity has on a spherically symmetric scalar held (i.e., 
Bray-Goetz (2014) [34], equation 10). The result is a scalar held whose spatial frequency 
becomes smaller as the distance from the origin increases, only to then drop off exponentially 
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Figure 1. Solution to the Klein-Gordon equation in a fixed spherically symmetric potential well 
based on the Milky Way Galaxy at t = 0, t = 10 Myr and t = 20 Myr. The pictures show the dark 
matter density (in white) in the xy plane. This solution, which one can see is rotating, has angular 
momentum. 

outside a given radius. However, for most of the region before the exponential drop off, the 
spatial frequency is quite constant, a property the solutions to equation 2.2 that we’ll be using 
have as well. While gravity will change these spatial frequencies, this is not a problem since 
we’ll be allowed to specify our spatial frequencies as we like later. Finally, we approximate 
the exponential drop off by simply declaring that / equals zero outside a given radius, chosen 
to be far outside the radius of the regular visible matter. 

Solutions to Eq. 2.2 can then be written as linear combinations of 

/ = Acos{ojt) • Yn{9,<t)) ■ r” • fuj,n{r) (2.3) 

where 

C,n{r) + = [m^ - UJ^)fuj,n (2.4) 

and Yn{9, 4>) is an nth degree spherical harmonic. Note that we require „(0) =0 but have 
not specified an overall normalization. Naturally, to get a complete basis of solutions we also 
need to include solutions like the one above but where cos{ujt) is replaced by sin(a;f). 

The form of the Minkowski solution on which our simulations are based is 

f = Aq cos{u:ot)f^gfi{r) + A 2 cos(a;2t - 2(/)) sin^(6')r^/a;2,2(r) (2.5) 

where /a;o,o(^) and satisfy Eq. 2.4. We note these solutions fall into the form of 

Eq. 2.3 since both cos(2(/>) sin^(0) and sin(2i?!)) sin^(0) are second degree spherical harmonics. 
Hence, the above Minkowski spacetime solution is the sum of a spherically symmetric solution 
(degree n = 0) and two degree two solutions. Note that these special solutions all have 180° 
rotational symmetry with densities which rotate rigidly. Note that this rigidity is not an 
assumption, but a consequence of approximating away the other spherical harmonic terms, 
assumed to be relatively small. As we said before, we conjecture that two armed spiral 
patterns result when terms similar to these dominate. 

Also, computations suggest that the approximation in Eq. 2.5 is a reasonable qualitative 
approximation to a solution in the fixed spherically symmetric potential case. For example, 
the pictures in Fig. 1 show the scalar field dark matter densities in the xy plane at t = 0, 
t = 10 million years, and t = 20 million years in a fixed spherically symmetric potential 
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well based on the Milky Way Galaxy. Like the Minkowski solution in Eq. 2.5, the solution 
depicted in Fig. 1 is also a sum of a spherically symmetric solution and two degree two 
solutions, but in a spherically symmetric potential. Notice how the solutions appear to have 
compact support, although actually the solution has simply decayed very rapidly outside a 
finite radius. Otherwise, the interference pattern is very similar to what we see with the 
solutions in Eq. 2.5, with the most obvious difference being in the last ring or two, which are 
a bit stretched out in the radial direction compared to our model. Since all we ultimately 
care about are the qualitative characteristics of the gravitational potential for the dark matter 
that we ultimately produce, this minor difference is something we are willing to tolerate to 
get simpler equations to solve. 

The reason the dark matter density is rotating in Fig. 1 and in our Minkowski model is 
most easily seen by making the substitution 


a = 



t 


into our expression for / in Eq. 2.5 to get 

/ = Aocos{ujot)fuiofi{r) + A2 cos{ujot- 2a) sin^(6')r^/^2,2(?')- 


( 2 . 6 ) 


(2.7) 


Note that to the extent that a stays fixed in time, then the above solution does not rotate 
and gives a fixed interference pattern since both terms are oscillating in time with the same 
frequency ujq. Hence, we get a wave dark matter interference pattern which is rotating 
according to the formula 

<>0 = (2.8) 

with period 

47r 

Tdm =-, (2.9) 

002 — OOq 

which may be related to the pattern period of the resulting barred spiral patterns in the 
regular matter. 

Plugging the wave function / defined in Eq. 2.7 into equation 2.1 defines both the wave 
dark matter density ^dm in the galaxy and the corresponding galactic potential x,y, z) 
due to the wave dark matter. Calculations using spherical harmonics as described in [27] then 
gives us 

~ Uo{r) + U2{r){2>z^ — r^) + Ui{r){2>hz‘^ — + Sr'^) 

+bi 2 (r)(r^ cos(2a) sin^(0)) (2-10) 

where 


56 

Uo{r) = Alf^,^,{rf + 

U2{r) = 

Ui{r) = 

U2{r) = 2AQA2fujafl{r)fuj2,2{r)- 
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Similarly, the rotating ellipsoidal galactic potential is then given by solving ode’s: 


V 

dTT/io 


Wo{r) + W2{r){3z‘^ - r^) + Wi{r){35z^ - + 3r^) 

+W 2 (r)(r^ cos(2a) sin^(0)), 


( 2 . 12 ) 


where for 0 < r < rmax, 

Kir)+‘^^'^^W:,ir) = Un{r), with <(0) = 0, lK(w) = 
and for r > rmax, 

Wn{r) = Wnirmax) • 

Note that the above conditions are meant to apply to W 2 {r) as well. 
Our input parameters for computing everything are: 


'^max 

2n + l 


(2.13) 


(2.14) 


• ^ 0 ) to be chosen to match the dark matter density in a given galaxy, 

• f'max, the radius of the wave dark matter density, 

• ^0 (typically = 1) and A 2 which control the relative proportion of the spherically 
symmetric and the second degree spherical harmonic components of the wave dark 
matter. 


• Wo and UJ 2 , the time frequencies of the spherically symmetric and the second degree 
spherical harmonic components of the wave dark matter, 

• m, a fundamental constant of nature in the Klein-Gordon equation, or equivalently the 
mass of the dark matter boson, yet to be determined precisely. 

From Eq. 2.4, we see that the spatial wavelengths of the spherically symmetric and the 
second degree spherical harmonic components of the wave dark matter are 


Afc = 


27r 


wf — 


(2.15) 


for k = 0,2. Then using the above equation and Eq. 2.9, we can solve for wq, W 2 , and m in 
terms of Aq, A 2 , and Tom to get 


Wo = 


TT 

2 

TT 


W2 = 


Ai 


Ai 




Ag 


= 




Tdm — 
Tdm + 

(^) 


27r 

Tdm 

2tx 


Tdm 
2 


Ai 


Ai 




Ai 


Ai 


Hence, if we prefer we may choose the input parameters to be 


(2.16) 

(2.17) 

(2.18) 


* /^O; ^maxi ^0) ^ 2 , Ao, A 2 , and Tdm, 
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Figure 2. Spiral Galaxy Simulation in [27]: The dark matter density times in the xy plane, 
the potential function in the xy plane, the level sets of the potential function in the xy plane, and the 
rotation curves. 


where Tqm is the time it takes the ellipsoidal galactic potential to rotate around once. These 
input parameters are more convenient in many instances and are probably more closely related 
to observable properties of galaxies. 

Bray (2010), in a series of test particles simulations, use the triaxial potential of Eq. 
2.12 to see the imprints that the wave dark matter, modeled in this way, produce on a disk of 
particles. Fig. 2, originally from that paper, shows how interference patterns in the wave dark 
matter may result in a rotating ellipsoidal galactic potential with plausible rotation curves. 


3 Comparison with a NFW triaxial halo 

As mentioned in the Introduction, there are in the literature some works dealing with the 
effects that a DM halo, with triaxial shape, has on the disk of its galaxy. For that reason 
here we also use a NFW halo to compare our results with the already existing studies. 

A NFW halo results from CDM simulations and is characterized by a density profile 
that rise steeply towards the center [35] 
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PNFwir) 


Ps 


(3.1) 


(r/rs)(l + r/rsY ’ 

with Ps and Vg as fitting parameters. 

Given a spherically symmetric density distribution / 3 (r), Chandrasekhar (1969) [36] gives 
a detailed procedure to generalize it into a triaxial distribution. Knowing p(r), the gravita¬ 
tional potential due to the generalized halo with triaxial shape is given by 


$ = 


fG p(i\a,20‘‘i 



(3.2) 


where ai, 02 and 03 are the semi axes of the triaxial distribution, and the function m^(tt) 
(not to be confused with the scalar field mass m) is defined as 




(3.3) 


By putting pnfw inside Eq. 3.2, the gravitational potential at an inner point of a NEW 
triaxial halo is given by 


^NFw{x,y,z) = iFGabcporg 

r°° m{u)/rs du 

^ Jo 1 -k m{u)/rs + u){b‘^ + u){c^ + u)' 

Finally we put this triaxial halo to rotate as a rigid body. 

We have two halo models with triaxial shapes accounted by very different approaches. 
It is also important to mention that although both kind of halos are triaxial, the dark matter 
nature proposed by each model (CDM and WDM) has implications on the density profiles of 
the halos, mainly on the inner regions. A NEW halo has a cuspy behavior whereas a WDM 
halo has a core behavior. So although both allow triaxiality, this will differentiate the models 
by using observational data, for example rotation curves. 

4 The code 

The simulations shown in this work were carried out using the latest version of the ZeusMP 
code [37] to model the baryonic component of a galaxy as a fluid which, as most of the material 
in an astrophysical system, has zero viscosity, very high compressibility, and low densities. 

By contrast, in [27] the stars, gas, and dust were modeled using test particles. One 
big difference with a fluid model is the requirement that there is a flow velocity vector u at 
each point which makes velocity a function of position, which is not the case with n-body 
simulations. However, in actual spiral galaxies, stars, gas, and dust do have very similar 
velocities when they are close to each other, which favors modeling them as a fluid in this 
regard. Hence, we consider the simulations presented here to be an important improvement 
over the previous simulations, plus being fully three-dimensional. 

With this hydrodynamic code the dark matter is implemented as an external gravita¬ 
tional potential in the gas momentum equation. This means that besides being affected by 
thermal pressure gradients, a fluid element in our simulations will feel the external influence 
of the gravitational potential in Eq. 2.12 or Eq. 3.4. 
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For each hydrodynamic simulation we set up the initial condition of a three-dimensional 
gas distribution as an isolated, cold, Miyamoto-Nagai disk with a density profile 

y^M aR^ -I- (a -|- 2>y/+ b‘^){a + V 

~ (R2 + (a + ^42^)2)5/2(^2 + 52)3/2 ' 

where M is the mass of the galaxy disk, a and b are the radial and vertical scale-length, 
respectively. 

We model the fluid as an ideal gas with internal energy density constant trough all the 
disk at the initial time. Fitting parameters to rotation curves data for LSB galaxies prefer 
thin disks, hence we chose a small value for the internal energy density ( 0.006J/m^), which 
means that the system is supported mainly by rotation, i.e., for the initial condition the gas 
component is distributed in the velocity space in such a way that it is supported by rotation 
against the background gravitational force due to the dark matter halo. In the next section 
we will see that the gas distribution does not undergoes a significant evolution other than a 
response to the triaxiality of the halo, with spiral patterns correlated to the semi axes of the 
distribution. 

5 Results 

First we test the potential models in Eqs. 2.12 and 3.4 against rotation curves data. As the 
main goal is to compare the two dark halos, we pick data from low surface brightness galaxies 
(LSB). An important feature of LSB galaxies is that the mass-to-luminosity ratio is usually 
higher than that of a normal spiral galaxy due to the low luminosity, and the dark matter 
fraction is much higher, which allows to model these systems as a unique component, a dark 
matter halo, neglecting the gravitational influence of the luminous matter. 

We set the parameters of the two models to fit the rotation curve data of HI in the LSB 
galaxy F568-3. This allow us to stablish the initial values in the simulations but we do not 
intend to model a particular galaxy. As we can see in Fig. 3, the rotation curve profiles are 
very different, specially in the inner parts, this reveals the cuspy nature of the NFW halo and 
the core nature of the WDM halo. 

With this initial fit we set the parameters of the two models (Eqs. 2.12 and 3.4) and 
evolve the hydrodynamics equations of a gas disk embedded in each triaxial halo. 

Figure 4 shows the results of our simulations. The luminous matter responds to the 
non-axisymmetric background, and although the triaxial shape in the halos is not achieved 
in the same way, is able to trigger the formation of a spiral pattern in both cases. 

In the particular case of WDM it is an important outcome, that triaxial halos arise 
naturally, a feature at galactic scales that WDM shares with CDM. This links the model in 
a natural way with the presence of spiral and barred patterns. 

Although the imprints of both halos on the gas are qualitatively the same, is worth 
trying to differentiate the underlying dark matter distribution, even if both are triaxial. 

To this purpose, since having the hydrodynamic simulations allow us to measure directly 
the rotational velocity of the gas, we outline the rotation curves over perpendicular directions 
to have an appreciation of the rotation of the gas and its possible dependence with the 
direction. 

Figure 5 shows the rotation curves of the gas embedded in the WDM (or SFDM) and 
NFW halos We measure the curves at a time in the simulation where non circular motions are 
evidently an important part of the gas kinematics. Blue and red lines indicate the directions 
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Figure 3. Initial rotation curve for the two triaxial halos adjusted to the data of the LSB galaxy 
F568-3. 



Figure 4. Face-on and edge-on projections of the gas distributions embedded in a WDM (or SFDM) 
triaxial halo and a NFW triaxial halo. In both cases the rotating triaxial potentials trigger the 
formation of spiral patterns. 


in which we outline these curves First we notice that the rotation curves in Figure 5 (bottom 
panel) change with the measurement direction, this changes are more notorious at large radii 
while in the central region seems to be smaller. 

On the bottom of those same plots we quantify this change by computing the difference 
between the two curves for each halo. 

In the case of the WDM halo we can see that AV is nearly constant in the inner Akpc, 
which means that although the curves are different, they rise with the same slope. 

For the NFW halo we see a similar behavior, AV is small and nearly constant within 
the inner 2kpc. Even more, for the central region the rotation curve along the two directions 
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Figure 5. (Top) Projected gas density, red and blue lines indicate the direction in which we measure 
the rotation curve. (Bottom) Rotation curves in two different directions for each halo. 


do not deviate significantly from the circnlar velocity (dotted line). 

The similar behavionr that the rotation cnrves have in the inner region of both halos is 
an indication that in spite the triaxiality of the backgronnd potential and the evident presence 
of non circnlar motions on the gas, observations conld still be able to distingnish the central 
distribntion of dark matter. 

The relevance of this resnlt for the case developed here and originally in Bray (2010) 
is that the natnral triaxiality of WDM halos preserve the inherent featnres of a scalar field 
confignration, a density core that becomes important when trying to explain the observed 
dark matter distribntion in dwarf galaxies. The spiral imprints that the triaxial halos have 
on the Inminons matter do not prevent the models of being tested by observations. 
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Figure 6. Face-on projection of the gas density embedded in a smaller, less massive, and less rounded 
halo whose parameters are adjusted to model the LSB galaxy F583-4. The gas develops open spiral 
arms and a bar-like structure. 

Finally we also adjust the model to describe the properties of a smaller and less massive 
LSB galaxy, the parameters in eq. 2.12 allow us to variate the triaxiality of the dark halo. In 
Figure 6 (top panel) we can see that the imprints of this less rounded halo on the baryonic 
matter is to stretch the gas distribution, which develops spiral arms with a greater pitch 
angle (more open arms) compared with the previous simulation. Also the disk central region, 
initially axisymmetric, develops a bar-like distribution. 

6 Discussion and Conclusions 

A large number of observations in LSB and dSph galaxies are better described by dark matter 
halos with cored density profiles, opposite to the predictions of N-body simulations with 
collisionless dark matter. A vast number of numerical studies have shown that this discrepancy 
has a plausible solution in the baryonic feedback which efficiently injects energy to the dark 
matter particles and therefore flattens the halo center. However the energy requirements 
and efficiencies for this mechanism to transform cusps into cores are not typical and mostly 
works at high redshifts [38, 39]. Also, simulations of dwarf galaxies show that even when 
the feedback produces a core in the center of the halo, a cusp could regrow by the infall of 
substructures to the galaxy [40]. 

Because the baryonic feedback solution seems to work only under specific circumstances 
and do not explains the presence of a core at all galactic scales, it is important to explore 
alternative DM models that predict cored density profiles. 

As in previous studies, here we explore the WDM model but generalizing the density 
profiles to rotating triaxial halos as this would be a more natural scenario. This is an impor- 
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tant outcome for the model, that in general the halos have triaxial shape, equivalent to the 
halos predicted by CDM simulations. 

The fact that the gravitational potential is no longer spherically symmetric or axisym- 
metric has important effects on the baryonic matter of the galaxy. It is not surprising to see 
the development of a spiral pattern in the galactic plane, which seems to be a usual outcome 
for a triaxial potential as shown in our comparison with a triaxial NFW halo. 

This is important because when trying to obtain information about the dark matter 
distribution in a galaxy, the mere presence of spiral arms or a bar may indicate that the 
baryonic matter dominates the central region and therefore observations like rotation curves 
may not tell us precisely what the DM distribution is at the halo center. But here we show 
that spiral arms and bars can be developed in a DM dominated galaxy with a central density 
core without supposing its origin on mechanisms intrinsic to the baryonic matter. These non- 
axisymmetric features are well explained and triggered by the interaction of the disk with the 
WDM triaxial halo, a triaxiality that arise naturally within the model. 
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